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Abstract 

Aerospace  has  routinely  processed  the  Oshornc  Time  Transfer  Receiver  (TTR)  data 
for  the  purpose  of  monitoring  the  performance  of  ground  and  GPS  atomic  clocks  in 
near  real-time  with  on-line  residual  displays  and  characterizing  clock  stability  with 
Allan  Variance  calculations.  Recently,  Aerospace  added  the  ability  to  estimate  the 
TTR’s  location  by  differentially  correcting  the  TTR’s  location  in  the  WGS84  reference 
system.  We  exercised  this  new  feature  on  a  set  of  TTR  clock  phase  data  and  obtained 
sub-meter  accurate  station  location  estimates  of  the  TTR  at  the  Aerospace  Electronic 
Research  Lab  (ERL). 

BACKGROUND 

The  Osborne  Time  Transfer  Receiver’s  (TTR)  primary  function  is  to  provide  a  means  to  monitor  a 
local  laboratory  frequency  standard  by  comparing  this  standard  to  GPS  system  time  using  the  GPS 
Navstar  satellites.  In  order  to  accomplish  this  task,  it  is  necessary  to  set  up  irutially,  the  local  receiver 
coordinates  into  the  receiver’s  memory.  This  initialization  is  assisted  by  using  the  built-in  positioning 
capabilities  of  the  receiver  whereby  real-time  location  estimates  are  obtained  during  user  specified 
tracks.  The  TTR  is  a  single  channel  Clear  Acquisition  (CA)  receiver,  thus  the  real  time  estimates  of 
location  are  obtained  sequentially  in  time  from  each  satellite.  This  operation  is  in  contrast  to  position 
receivers  designed  to  receive  signals  from  four  or  more  satelHtes  simultaneously  using  the  Precise  (P) 
Code,  which  is  an  order  of  magnitude  more  accurate  than  the  CA  code.  The  lack  of  simultaneity  of  the 
measurements  further  degrades  the  position  estimates  from  the  TTR.  Therefore  the  initially  entered 
station  coordinates  are  only  first  estimates  and  subject  to  error. 

These  deficiencies  are  unessential  to  the  receiver’s  operation  as  a  time  transfer  device,  and  consequently 
the  TTR  has  been  most  useful  and  helpful  in  quickly  detecting  atomic  clock  anomalies  not  only  in  the 
laboratory  standard,  but  also,  satellite  atomic  clocks  themselves.  This  latter  ability  is  due  to  the  high 
precision  and  stability  of  the  laboratory  clock. 

The  initially  entered  receiver  locations  were  obtained  by  observing  the  TTR’s  solution  on  a  daily  basis 
for  several  days.  Since  no  consistent  daily  solution  could  be  obtained,  the  locations  were  finally  selected 
and  entered  by  a  trial  and  error  method.  The  receiver’s  initial  position  errors  are  clearly  manifested 
in  the  clock  phase  difference  measurements  as  pronounced  systematic  diurnal  signatures  in  some  the 
satellites,  (Satellites  ID  PRN  9  and  PRN  3)  as  seen  in  figure  1.  The  long  non-diurnal  systematic  effects 
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observed  for  each  satellite  seen  in  this  figure  is  the  effect  of  clock  steering  (confirmed  by  independent 
data  from  USNO)  performed  by  the  GPS  control  segment.  The  clock  steering  is  necessary  to  keep 
GPS  system  time  to  within  a  microsecond  of  USNO  time. 

The  program  which  was  used  to  monitor  atomic  clock  performance  and  characterize  clock  stability 
by  computing  the  Allan  Variance  was  modified  with  additional  code  for  an  iterative  least  squares 
procedure.  The  initial  estimate  is  refined  by  inputting  the  time  difference  measurements  from  the 
receiver  as  data  into  this  least  squares  program  and  thus  derive  a  better  estimate  for  the  station  coor¬ 
dinates.  Since  the  original  station  entries  are  referenced  to  the  WGS84  geodetic  coordinates  as  given 
by  the  TTR,  and  the  broadcasted  ephemerides  are  also  based  on  the  WGS84  system,  the  resulting 
least  squares  corrections  will  be  referenced  to  the  WGS84  coordinates. 


ANALYSIS 

We  note  that  the  Ar,  the  range  difference  (i.e.time  difference  measurement  error  times  the  speed  of 
light)  is  merely  the  negative  sum  of  the  projections  of  the  South,  East  and  height  components  of  the 
station  error  into  the  line  of  sight  from  the  station  to  the  satellite.  A  derivation  of  this  is  given  in  the 
Appendix  section. 


Ar  =  ~  cos  AZ  cos  e/(  South  error)  —  s\nAZ  cose/(  East  error)  —  sin  e/(  Height  error)  (1) 

where  (e/)  and  {AZ)  are  the  observe  elevation  and  azimuth  angles.  Azimuth  is  defined  as  measured 
from  North  round  by  East,  or: 


Ar  =  —  cos  AZ  cos  el{RA<j>)  —  sm  AZ  cos  el{Rcos  (f>AX)  ~  sine/(A/i)  (2) 

where  Ah,  A0  and  AA  are  the  respective  height  correction  above  a  reference  geoid,  and  corrections 
in  latitude  and  longitude  of  the  TTR.  The  system  of  normal  equations  necessary  for  a  least  squares 
solution  is  given  by  equation  (3),  where  the  left  hand  side  contains  the  data,  and  the  right  hand  side 
consists  of  the  station’s  South,  East  and  height  errors  to  be  determined  from  the  data.  Each  row 
of  the  matrix  of  normal  equations  ||A||  is  formed  from  the  coefficients  of  equation  (2).  A  vector  of 
measurements,  !Ar|,  is  formed  from  the  data  set,  Ar.  The  solution  proceeds  in  an  iterative  two  step 
procedure.  The  data  are  first  fitted  to  a  polynomial  up  to  sixth  degree  in  time  to  remove  the  clock 
phase  wander  and  other  systematic  effects.  The  remaining  difference,  Ar,  is  used  in  a  standard  least 
squares  solution  for  the  station  off-sets  in  the  South,  East,  and  height  directions,  as  follows: 


^  (South  error)  \ 

(East  error)  =  [||A|f  ||A||]  ||A|f  |A7^  (3) 

(Height  error)  j 

As  many  as  15  iterations  are  necessary  to  converge  the  final  solution.  Each  iteration  includes  a 
reestimation  of  the  polynomial  to  separate  out  the  clock  wander  and  steering  effects  from  the  station 
location  estimation. 
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DATA  DESCRIPTION  AND  RESULTS 


A  span  of  clock  phase  difference  data  between  2/8/90  to  d/lO/OO  from  the  TTR5  was  selected  for 
the  receiver  coordinate  adjustment.  This  data  also  contained  the  elevation  and  azimuth  angles  as 
computed  by  the  TTE,5,  which  are  needed  in  equation  (2).  The  iterative  least  squares  correction  for 
the  coordinates  is  obtain  by  computing  the  weighted  correction  from  each  NAVSTAR.  At  the  end  of 
each  iteration,  the  sum  of  the  average  data  residual  (standard  deviation  or  s.d.)  from  each  NAVSTAR 
is  computed,  and  the  iteration  ceases  when  this  sum  reaches  a  minimum.  The  result  and  solution  from 
this  procedure  are  summarized  in  Table  1. 

As  seen  in  this  table,  the  final  post-fit  data  s.d.  from  each  NAVSTAR  are  substantially  reduced  from 
the  initial  pre-fit  data  s.d.  The  columns  labelled  “corrections”  are  the  final  corrections  determined 
from  data  from  each  NAVSTAR,  and  show  considerable  scatter  from  NAVSTAR  to  NAVS'fAR  (e.g. 
-11.726  meters  of  West  error  for  NAVSTAR  9).  The  individual  s.d.  for  each  NAVSTAR  for  each  coor¬ 
dinate  component  are  listed  under  the  column  labelled  “sig”.  The  final  solution  for  the  South,  East, 
height  correction  components  and  the  associated  composite  s.d  from  all  the  NAVSTARs  are  listed 
under  the  columns  labelled  “weighted  mean”  and  “composite  sigma”  respectively.  The  corrections  are 
7.366m  South,  5.612m  West  and  16.1 2m  low. 

These  formal  s.d.  need  to  be  in  ultiplied  by  a  factor  of  about  6  to  7  corresponding  to  the  postfit  average 
data  residual  as  described  by  the  data  s.d.  The  covariance  implicit  in  equation  (3)  assumes  a  data 
sigma  of  one  meter,  whereas  it  is  more  appropriate  to  use  the  actual  data  sigma  of  about  7  meters. 
Consequently,  the  more  realistic  uncertainties  are  0.60  meters,  0.39  meters  and  0.36  meters  in  the 
respective  coordinates.  Finally,  when  the  linear  corrections  are  converted  to  arc  second  measure,  the 
corrections  are  0.2382”  South  in  latitude  and  0.2187”  West  in  longitude. 

When  these  corrections  are  used  in  equation  (2)  to  compensate  the  data  for  station  error,  a  marked 
improvement  is  seen  when  Figure  2  is  compared  to  Figure  1.  The  most  noticeable  effect  is  the  reduc¬ 
tion  in  the  diurnal  phase  errors  for  NAVSTAR  PRN  9  and  3  as  well  as  moderate  error  reduction  for 
the  other  NAVSTAR  FRNs. 

For  further  confirmation  of  the  validity  of  these  results,  we  performed  a  “blind  test”  verification  by 
using  the  corrective  factors  derived  from  data  between  2/8/90  to  4/19/90  into  a  data  span  not  used 
in  the  solution.  We  selected  a  earlier  set  between  10/3/89  to  11/12/89.  The  result  of  this  test  is 
shown  in  Figures  3  and  4.  Figure  3  shows  the  errors  when  the  station  correction  is  left  out.  As  can  be 
seen,  sizable  errors  for  PRNs9  and  3  are  there  and  of  the  same  magnitude  as  seen  in  Figure  1,  where 
the  station  corrections  are  also  left  out.  In  P’igure  4,  where  the  corrections  are  applied,  significant 
reductions  in  the  diurnal  errors  especially  for  PRNs  9  and  3  are  achieved.  As  an  additional  note,  the 
data  selected  were  time  differenced  measurements  between  the  Aerospace  laboratory  standard  and 
each  of  the  NAVSTAR  clocks.  Hence,  as  can  been  seen  in  Figures  3, and  4,  the  individual  satellite 
clock  behavior  is  exhibited.  A  summary  of  each  of  the  clock’s  frequency  labeled  as  FRQ  and  aging 
labeled  as  DFQ  between  each  satellite  and  the  laboratory  is  shown  in  the  upper  right  hand  corner  of 
the  figure.  The  units  for  frequency  difference  between  laboratory  standard  and  satellite  clock  are  in 
picoseconds/second  and  the  aging  factor  are  in  nanoseconds/days'^. 
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SUMMARY  AND  CONCLUSION 


A  few  meters  of  TTR  location  input  error  will  result  in  noticeable  effects  in  the  data.  Significant  im¬ 
provements  resulted  when  station  corrections  of  0.2382”  South  in  latitude,  0.2187”  West  in  longitude 
and  -16.12  meters  in  receiver  height  are  applied.  The  results  here  show  that  the  TTR  with  CA  coded 
data  provides  sufficient  precision  for  the  receiver  to  act  as  a  WGS84  position  locator  with  sub-meter 
accuracy  as  well  as  a  time-transfer  device. 

APPENDIX 

For  the  purpose  of  differentially  correcting  the  station  coordinates,  we  formulate  the  rectangular  station 
components  for  a  spherically  shaped  earth, 

X  =  h)  cos  ^  cos  X 
Y  —  (R  +  h)  cos  4>  sin  A 
Z  =  {RYh)  sm<j} 

where  R,  h,  4>  and  A  are  the  respective  radius,  height  above  a  reference  geoid,  geocentric  latitude  and 
longitude  of  the  TTR. 

We  define  a  set  of  unit  vectors  L,  A,  D  as: 

Lx  =  cos  4>  cos  A  Ai;  =  —  sin  A  Dx  =  —  cos  A  sin  <f) 

Xy  =  cos  </)  sin  A  Ay  =  cos  A  Dy  =  —  sin  A  sin  (1) 

Xj  =  sin  (f)  Ax  =  0  Dz  —  cos  (p 

so  that  errors  in  the  rectangular  components  are  related  to  errors  in  the  spherical  components  through 
the  matrix  equation  as: 


f  AX  \ 

(  Lx 

Ax 

Dx 

] 

( 

Ah 

\ 

AY 

= 

Ly 

Ay 

Dy 

R  cos  d>AA 

\  ) 

^  Lz 

Dz 

\ 

RAP 

/ 

In  order  to  obtain  the  errors  AXh,  AV)j,  and  AZ^  in  a  local  horizon  tangent  plane,  with  the  normal 
to  the  plane  in  the  local  zenith  direction,  we  need  to  rotate  AX,  AT,  AZ  by  a  rotation  matrix: 


1  AXh  \ 

(  -Dx 

~Dy 

-Dz 

AYh 
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\  AZh  ^ 
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Substituting  AX,  AT,  AZ  from  equation  (2)  into  (3),  we  obtain: 


1  ] 

^  —Dx  —Dy  —Dz 

^  Lx  Ax  Dx  \ 

'  Ah  \ 

AYh 

— 

-^y  -^Z 

Ly  Ay  Dy  1 

R  cos  pAX 

\  j 

\  / 

y  Lz  Az  Dz  J 

y  RAp  y 

which  reduces  simply  to: 


/  AXh  \ 

/  0  0  -1  \  / 

AYh 

= 

0  10 

y  AZh  y 

o 

o 

Ah  \ 
i2cos(^AA 
RA(t>  ) 


(3) 


(4) 


(5) 
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The  topocentric  rectangular  measurement  error  vector,  Ap,  with  components,  An,  An,  Aw,  is  the 
negative  of  the  station  error  vector  AR,  (assuming  no  errors  for  the  NAVSTARs’  position). 

A^Au,  An,  Aw)  =  -AR{AXh,  AYh,  AZh)  (6) 

A  triad  of  unit  vectors  in  the  topocentric  system  is  defined  for  the  observed  elevation  (e/)  and  azimuth 
{AZ)  angles  as  follows: 


so  that: 


=  —  cos  el  cos  AZ  A^h  =  sin  AZ 
Lyh,  =  cos  el  sin  AZ  Ayh  =  cos  AZ 

L^h  =  sin  el  A^h  =  0 


Dxh  ~  cos  AZ  sin  el 
Dyh  =  -  sin  AZ  sin  el 
Dzh  —  cos  el 


^  Aw 

^  Lxh  Dxh  ^ 

(  At  ^ 

Aw 

= 

Lyh  ^yh  Dyh 

r  cos  elAAZ 

^  / 

^  Lzh  Azh  Dzh  j 

^  rAel  j 

The  inverse  is  then: 


^  Ar 
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1  Lxh 

Lyh 

Lzh 

1 

Aw 

\ 

r  cos  elAAZ 
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(8) 


Noting  the  relationship  between  Ap  and  AR  in  equation  (6),  and  substituting  (5)  into  (8)  we  get  for 
first  component,  the  range  expression,  Ar  as: 


Ar  —  -  cos  AZ  cos  el(RA(f>)  -  sin  AZ  cos  el(R  cos  4>AX)  ~  sine/(A/i)  (9) 


or: 


Ar  =  —  cos  AZ  cose/(  South  error)  —  sin  AZcose/(  East  error)  -  sine/(  Height  error)  (10) 
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AEROSPACE  ELECTRONIC  LAB  TTR 
LEAST  SQUARES  STATION  ADJUSTMENT 
(meters) 


ITERATIVE  SOLUTIONS  FOR  EACH  NAVSTAR  PRN 


stnd  deviation 


pre-fit 

post-fit 

final  South 

final  East 

final  Height 

residual 

residual 

correction 

correction 

correction 

PRN 

correction 

sig 

correction 

sig 

correcton 

sig 

6 

6.809 

6.415 

-2.776 

1.115 

-  1.013 

0.971 

1.667 

0.568 

9 

9.298 

7.350 

-0.215 

0.116 

-11.726 

0.801 

-0.624 

0.083 

11 

7.146 

6.557 

-6.023 

0.430 

5.883 

0.410 

1.002 

0.150 

12 

6.044 

5.888 

1.054 

0.203 

-2.926 

0.209 

-0.768 

0.108 

13 

4.640 

3.691 

-2.125 

0.331 

-0.064 

0.081 

1.285 

0.206 

3 

7.690 

6.570 

1.857 

0.179 

0.234 

0.081 

-0.661 

0.102 

WEIGHTED  SOLUTION 


South  Correction 
weighted  composite 
mean  sig 

-7.366  0.600 


East  Correction 
weighted  composite 
mean  sig 

-5.612  0.387 


Height  Correction 
weighted  composite 
mean  sig 

-16.122  0.360 


AEROSPACE  ELECTRONIC  LAB  TTR 
STATION  LOCATION 


Latitude 

Longitude 

Height 

Old  value: 

-|-33deg 

54min 

54.5778sec 

241 deg 

37min 

13.5410sec 

25.00  meters 

Correction: 

-0.2382sec 

-0.2187sec 

-16,12  meters 

New  value: 

+33deg 

54min 

54.3396sec 

241deg 

37min 

13.3223sec 

8.78  meters 

TABLE  1. 

Summary  of  solutions  and  estimates  of  TTR5 
station  location  in  WGS84  reference  system 
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Figure  1.  Clock  phase  residuals  between  ground  cesium  and  GPS  system  time 
via  each  NAVSTAR  PRN.  Same  systematic  effect  for  satellite  is  due  to  GPS  time 
steering.  Diurnal  “scatter”  is  due  to  station  location  error. 


Figure  2.  Clock  phase  residuals  between  ground  cesium  and  GPS  system  time 
via  each.  NAVSTAR  PRN.  Diurnal  “scatter”  considerably  reduced  from  Figure  1. 
Reduction  is  due  to  corrected  station  location. 
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Figure  3.  Clock  phase  residuals  between  ground  cesium  and  each  NAVSTAR 
PRN.  The  individual  satellite  clock  variations  are  apparent.  The  table  on  the  up¬ 
per  right  corner  of  figure  lists  the  frequency  offset  (FRQ)  in  units  of  picosec/sec. 
and  ageing  (DFQ)  in  units  of  nanosec/days^.  Diurnal  “scatter”  is  due  to  station 
location  error. 
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Figure  4.  Clock  phase  residuals  between  ground  cesium  and  each  NAVSTAR 
PRN.  The  table  on  the  upper  right  corner  of  figure  lists  the  frequency  offset  (FRQ) 
in  units  of  p;cosec/sec.  and  ageing  (DFQ)  in  units  of  nanosec/days^.  Diurnal  “scat¬ 
ter”  considerably  reduced  by  using  same  station  locations  values  as  used  in  Figure 
2.  and  listed  in  Table  1. 


